Parameterization of Dark-Energy Properties: a Principal-Component Approach 
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Considerable work has been devoted to the question of how best to parameterize the properties 
of dark energy, in particular its equation of state w. We argue that, in the absence of a 
compelling model for dark energy, the parameterizations of functions about which we have no 
prior knowledge, such as 'w{z), should be determined by the data rather than by our ingrained 
beliefs or familiar series expansions. We find the complete basis of orthonormal eigenfunctions 
in which the principal components (weights of 'w{z)) that are determined most accurately are 
separated from those determined most poorly. Furthermore, we show that keeping a few of the 
best-measured modes can be an effective way of obtaining information about w{z). 



The discovery of the accelerating universe P, Q has 
been followed by a lot of work concentrated on how best 
to parameterize dark energy and measure its properties 
(d; ^ IE Hi Q and references therein). The parameters 
of choice have been the energy density in this compo- 
nent, flx: and its equation-of-state ratio w{z) Q. While 
current data allow interesting constraints only if w is as- 
sumed constant, future surveys, such as those envisioned 
for the Supernova Acceleration Probe (SNAP) ^] or the 
Large-aperture Synoptic Survey Telescope (LSST) ^ 
may allow relaxing that assumption and measuring, for 
instance, the first derivative of w with redshift, w' = 
dw/dz\z=o. Despite these exciting prospects, it has be- 
come clear that accurate constraints on a general w{z) 
will remain elusive even for the most ambitious sur- 
veys ^3:, 4]; hence, a simple but general parameterization 
of w{z) seems necessary. 

Yet some issues remain obscure; in particular the ques- 
tion of what is actually being measured by a given sur- 
vey. In other words, what are the quantitative weights 
on measurements of w{z)l Consider, for example, a 
fiducial dark-energy model with the equation of state 
w{z) = — l-|-sin(7rz). Measurements of w' would indicate 
a positive quantity even though the redshift-averaged 
derivative of w is zero; this is because the largest sensi- 
tivity of a measurement is at the low-redshift end where 
w'{z) > 0. Using somewhat more general arguments, 
it has been argued 0, Q that w{z) is best determined 
around z ~ 0.3, and less accurately determined at much 
lower and much higher redshift because of the Hubble law 
and decreasing importance of dark energy respectively. 
We would like to better understand this and, moreover, 
find a natural basis of weights that represents the mea- 
surement. 

We consider the coordinate distance r(z) as the pri- 
mary observable of any survey (although this could eas- 
ily be generalized to number-count measurements, for 
example) and assume a flat universe with energy den- 
sity in matter flm = 0.3 and that in the dark component 
ilx — We divide the redshift range of the survey 

into N bins centered at redshifts Zi with corresponding 



widths Azi [i = 1, . . . , N). We assume dark energy is 
parameterized in terms of w{z), which we define to be 
constant in each redshift bin, with a value Wi in bin i. 
In the limit N oo, this allows for a completely general 
w{z). We use N — 50, which provides sufficient reso- 
lution yet doesn't require large computational time (we 
have explicitly checked that the results change little for 
N > 20). For piecewise-constant 'w{z), the energy den- 
sity of the dark component evolves as (for z in bin j, 
Zj - Azj/2 < z < Zj + Azj/2) 



px{z) ^ px{z ^ 0) 
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Then H\z) = [{1-VIm)px{z)/ px{OH^M{l + zf] 
and Hor{z) = J Hq/ H{z')dz' . Note that one could use 
the same methodology to reconstruct other cosmological 
functions, for example f{z) = Hq/H(z). 

In this paper, we only consider the function w(z). Al- 
though f{z) is easier to measure due to the fact that it 
is more directly related to the observable luminosity dis- 
tance, in order to understand the behavior of dark energy, 
one needs to measure its evolution with time and the 
scale factor — therefore, measure a derivative of pxjz) 
or f{z). Since w{z) is related algebraically to df/dz 
one effectively needs w{z). Therefore, despite the fact 
that measurements of f{z) or px{z) will be more accu- 
rate than those of w(z), constraints on the latter quantity 
will be crucial for understanding the nature of dark en- 
ergy. 

We use the Fisher matrix formalism to compute the co- 
variance matrix for the parameters Wi. We also include 
the "nuisance parameter" M., which is an overall offset 
in the magnitude-redshift diagram [l| and we then 
marginalize over these parameters. For the fiducial sur- 
vey, we assume 3000 type la supernovae (SNe) uniformly 
distributed in redshift between z = and z = 1.7. Any 
actual survey distribution will undoubtedly be somewhat 
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FIG. 1: The principal components of w{z) for Hm per- 
fectly known. The four best-determined and two worst- 
determined eigenvectors are shown and labeled for clar- 
ity. We have marginalized over the magnitude offset M. 

non-uniform in redshift, but we ignore this in order avoid 
seeing the effects due to the distribution itself. 

Principal components of the Fisher matrix. We start 
by computing the Fisher matrix F for our parameters Wi 
{i — I, . . . , N), having marginalized over A4 and assum- 
ing for a moment that ^Im has been determined inde- 
pendently. We find the familiar result that the accuracy 
in measurements of Wi , \/{F^^, rapidly increases with 
increasing i (that is, redshift). Now, it is a simple matter 
to find a basis in which the parameters are uncorrelated; 
this is achieved by simply diagonalizing the inverse co- 
variance matrix (which is in practice computed directly 
and here approximated with F). Therefore 

F = W^AW (2) 

where the matrix A is diagonal and rows of the decorrela- 
tion matrix W are the eigenvectors ei{z), which define a 
basis in which our parameters are uncorrelated . The 
original function can be expressed as 

JV 

■u;(z) = ^ Qfi ei(z) (3) 

4=1 

where Ci are the "principal components" . Since F is real 
and symmetric, W can be chosen to be orthogonal and 
of unit determinant, so that {ei{z)} form an orthonor- 
mal basis. Using the orthonormality condition, the coef- 
ficients ai can be computed as 

JV 

ai ='^w{za)e^{za). (4) 

a=l 

Diagonal elements of A, A^, are the eigenvalues which 
determine how well the parameters (in the new basis) can 



be measured; a{ai) = X~ . We have ordered the a's so 
that a{ai) < a{a2) < ■ ■ ■ < a{aN)- 

The principal components of F are shown in Fig. ^ 
For clarity, we only show 4 of the best-determined and 
2 of the worst-determined eigenvectors for w{z) for the 
fiducial model w(z) = —1. The best-determined com- 
ponent of w peaks at z w 0.2, enters with the coeffi- 
cient ai = —4.18, and can be measured to an accuracy 
of 2.5%. We find that the eigenvectors depend weakly 
on the fiducial cosmological model, while they depend 
somewhat more strongly on the set of cosmological pa- 
rameters; for example, marginalizing further over Q,m 
with a prior of 0.03, the best-determined eigenvector of 
w peaks at z w 0.1 while the others are mostly unaf- 
fected. In addition, allowing that Qm is not precisely 
known increases the uncertainty of how well the first few 
modes can be measured by 50-100%, while not affect- 
ing the higher modes much. Finally, note that, although 
the best-measured modes peak at relatively low redshifts, 
one still needs the full redshift range (z < 2) in order to 
measure those modes accurately. 

Note a few nice things about this decomposition (which 
has recently also been applied to the weak lensing re- 
construction of the radial density field 0|). First of 
all, the measured eigenvectors are orthonormal and or- 
dered by how accurately they can be measured. The M^^ 
best-determined eigenvector has precisely M — 1 nodes, 
which makes the interpretation quite natural: the first 
eigenvector corresponds to the "average of w{z)" , the 
second one to the "first derivative of w" , the third one 
to the "second derivative of w", etc. We argue that, 
in the absence of a compelling theoretical argument for 
any particular parameterization of 'w{z) except perhaps 
w{z) = —1, these modes are the natural basis in which 
the function w{z) should be considered for any given ex- 
periment. Furthermore, the spread in accuracies is larger 
than before: the best-determined modes are better deter- 
mined than the best-determined (original) components 
Wi, while the worst-determined modes are determined 
more poorly than the worst-determined Wi. Of course, 
the total constraint upon w{z) at any given redshift Za 

/ JV \i/2 
a{w{za))^[J2<J^a.)eUza)j (5) 

does not depend on the chosen basis. Consequently, by 
exploiting only a few of the best-determined modes we 
might hope to gain in accuracy while not biasing the 
result too much; we discuss this later on. 

A "sweet spot"? The best-determined eigenvector 
of w{z) peaks at z ~ 0.2, although other parameter 
choices (not including the magnitude offset A4, for ex- 
ample, or marginalizing over Qm) can change this to 
other values, even z — > 0. At first this may seem at odds 
with the result that w{z) is typically best determined 
around z ~ 0.3 However, the question of where 
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the "sweet spot" - the point of minimal uncertainty - 
of w{z) falls depends on the fmiction which is used to 
parameterize w{z). For example, using the parameter- 
ization u;(z) = wq + w' z leads to a sweet spot around 
z ~ 0.3. Fitting a low-order polynomial to the distance 
data may lead to one or two sweet spots in the fitted func- 
tion However, the equation-of-state ratio cannot 
truly be isolated at a given redshift; moreover, there is 
no compelling theoretical motivation for any particular 
parameterization. We therefore argue that the weights 
we compute are a better representation of the sensitivity 
of w{z) than the sweet spot in any particular parameter- 
ization. 

Testing the Constancy ofw{z). The only functional 
form for w{z) for which, we would argue, there is a strong 
theoretical bias is a constant w, in particular w{z) = — 1. 
It is therefore important to test whether the equation 
of state ratio w is constant or not. Non-constant w{z) 
would point toward a dynamical mechanism for late-time 
acceleration of the universe. To test constancy of w one 
could use actual r(z) data to compare for models of 
constant w{z) vs. those of varying Or one could 

simply measure w' [3ll3.ll4j. Here we seek a more general 
approach. 

Although the Mth eigenvector eAf(z) roughly corre- 
sponds to the {M — l)th derivative of wiz), it unfortu- 
nately cannot be used directly to test the constancy of 
w simply because the coefficient aM is not zero even for 
w{z) = const. Now, from Eq. |0J it follows that, for 
constant w(z) = Wc 

N 

di = Wc ^ei(za). (6) 

a=l 

where {di} correspond to the hypothetical Wc- One can 
then perform a simple t^st to determine whether w is 
constant: 

where we have chosen to keep only the first M eigen- 
modes, since the best-determined modes will contribute 
the most to the sum (the test is valid regardless of the 
value of M). Since the best-fitting Wc is to be deter- 
mined by finding the minimum there are M — 1 
degrees of freedom. For example, for a fiducial model 
w{z) = —1.0 -|- 0.3 keeping M = 4 modes rules out 
constant w at 98.7% CL. For comparison, the standard 
wq-w' test applied to this fiducial model rules out con- 
stant w at about 98.5% CL, since a{w') « 0.13. Al- 
though the two tests give comparable results, the former 
describes w with four parameters rather than two and is 
therefore more general than the w' test. In particular, 
the proposed test has a nice feature that the number of 
principal components that are used, M, can be chosen so 
as to maximize the strength of the desired constraint. 



Function Reconstruction? From Fig. ^ it is apparent 
that the well-determined eigenmodes: (1) are non-zero at 
low end of the redshift range and zero at the high end, 
and (2) do not oscillate much. The opposite is true for 
the poorly determined eigenvectors. This tells us that 
the accuracy in determining w{z) is best if this function 
is smooth and if we are trying to determine it at low to 
moderate redshift (z ^ 1). 

We have seen above that the eigenmodes of w{z) are 
roughly ordered by the absolute size of their coefficients 

- more noisy modes contribute less to w(z). One can 
then use the first M eigenvectors (those measured most 
accurately) in order to approximately reconstruct w{z): 

N M 

w{z) ^^aiei{z) Ki^aiei{z) (8) 

i=l j=l 

where Af < N. Obviously, in the case of M = TV we 
recover the original error bars, which are typically hope- 
lessly large (at least for N > few, which gives sufficiently 
high resolution in redshift). When M < N, two things 
happen: w{z) is reconstructed less accurately (so the bias 
increases) , but the error bars are smaller (so the variance 
decreases). Indeed, getting rid of the noisy eigenmodes 
corresponds to setting w{z) at high redshift end to zero, 
which will bias any w{z) that doesn't actually go to zero 
at that redshift. 

Fig. 121 illustrates these arguments (we choose = 20 
for definitiveness). To choose the optimal number of 
eigenmodes to be kept, M the statistically correct thing 
to do is to minimize risk [l5j, where 

risk = bias^ + variance (9) 

N N 

= ^ iw{z,) - w{z,)f + ^ a'iw{z,)), (10) 

i=l i=l 

the sums in Eqs. l|2Jl and |SJ| run from 1 to A/, and w is the 
fiducial value of w. The top panel of Fig.|21illustrates the 
minimization of risk. The bottom left panel of the Figure 
shows that the fiducial w{z) = 1/2(^1 -I- tanh(2 — 0.5)), 
which smoothly asymptotes to zero at z ^ 0.8, can be 
reconstructed without much bias; its risk is minimized by 
keeping Af = 4 terms. The same is not true for w{z) = 

— 1 + 0.3 z (bottom right; Af = 5 for minimum bias), 
which doesn't go to zero within our redshift range and is 
therefore reconstructed with a large bias. Furthermore, 
the reconstructed quantities that are linearly related to 
w(z), such as l+w{z), will also go to zero at high redshift 
a M < N (bottom right panel). The resulting bias in 
w{z) is therefore different for the two cases, and depends 
on whether w{z) or, say, 1 -I- w{z) is reconstructed. 

This analysis shows that, while reconstructing w{z) at 
z ^ 1 is somewhat promising, it does not seem possible 
to recover w{z) at z > 1 in a truly model-independent 
way. We have checked that even adding a large number 
of SNe at the high-redshift improves the reconstruction 
only marginally, leaving a significant bias at the high red- 
shift end. The fact that bias (and therefore risk) will be 
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FIG. 2: Reconstruction of w{z) by keeping only a fraction 
of eigenvectors so as to minimize risk. Top panel: illus- 
tration of the minimization of risk. Bottom left: optimal 
reconstruction (68% C.L. error bars shown) of fiducial 
w{z) (solid line) that goes to zero at high-redshift end. 
Bottom right: optimal reconstruction of w{z) that does 
not go to zero at high-redshift end. Also shown on this 
panel is optimal reconstruction of 1 -I- w{z) for the same 
w{z) model. 



difficult to estimate in practice exacerbates the problem. 
Note too that the proposed reconstructions of the equa- 
tion of state ratio are parametric, since they use 
a fitting function to fit the distance-redshift data. While 
planned surveys may provide accurate determinations of 
the fitting parameters, they cannot test well the validity 
of the survey parameterization. 

Conclusions. A number of models that roughly ex- 
plain the observable consequences of dark energy were 
proposed in recent years, some of them starting from fun- 
damental physics and others being purely phenomcnolog- 
ical. However, the origin and nature of dark energy re- 
main unknown, and it is safe to say that none of the mod- 
els should be taken too seriously at this point. Therefore, 
it seems wise to approach the question of constraining the 
properties of dark energy empirically, with as few prior 
assumptions as possible. 

Here we propose that rather than using various pa- 
rameterizations proposed in the literature to describe 
dark energy (which typically parameterize the equation 
of state w{z) using series expansions etc.), one can simply 
let data decide which weights of these functions are mea- 
sured best, and which ones are measured most poorly. 
These weights are the natural basis that parameterizes 
the measurements of any particular survey. 

For definitiveness, we assumed a cosmological distance- 
redshift survey containing 3000 SNe uniformly dis- 
tributed in redshift. We computed the weights of w{z) 



and showed that accurately measured modes (weights) 
are rather smooth and go to zero at higher redshifts, 
while the opposite is true for poorly measured modes. 
The previously-considered "sweet spot" in the sensitivity 
of w{z) is largely a function of the choice of parameteri- 
zation, and the shape of the first principal component is 
a better indicator of the redshift (s) at which w is being 
measured. 

With the proposed parameterization, the test of 
whether w{z) is constant is straightforward and intu- 
itive. Although the reconstruction of w{z) is straight- 
forward to implement, the reconstructed w{z) is noisy, 
and it is advantageous to keep only the best-measured 
modes in order to decrease the statistical reconstruction 
error. This introduces a systematic bias at z 1 for 
most models, roughly independently of the redshift cov- 
erage of SNe. We therefore conclude that while model- 
independent statements about w{z) at z < 1 may be 
feasible, those at z > 1 will be unreliable. 

In our opinion, the greatest advantage of this approach 
is simply having an intuitive and quantitative answer as 
to what is actually being measured by a given survey. 
The next logical step is to apply this method to other 
cosmological tests. We will address this in a future pub- 
lication. 
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